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ABSTRACT 

We present a Spitzer mid-infrared survey of 42 Fanaroff-Riley class II ra- 
dio galaxies and quasars from the 3CRR catalog at redshift z < 1. All of the 
quasars and 45 ± 12% of the narrow-line radio galaxies have a mid-IR luminosity 
of uL u (15^m) > 8 x 10 43 erg s _1 , indicating strong thermal emission from hot 
dust in the active galactic nucleus. Our results demonstrate the power of Spitzer 
to unveil dust-obscured quasars. The ratio of mid-IR luminous narrow-line radio 
galaxies to quasars indicates a mean dust covering fraction of 0.56 ±0.15, assum- 
ing relatively isotropic emission. We analyze Spitzer spectra of the 14 mid-IR 
luminous narrow-line radio galaxies thought to host hidden quasar nuclei. Dust 
temperatures of 210-660 K are estimated from single-temperature blackbody fits 
to the low and high-frequency ends of the mid-IR bump. Most of the mid-IR 
luminous radio galaxies have a 9.7 fim silicate absorption trough with optical 
depth < 0.2, attributed to dust in a molecular torus. Forbidden emission lines 
from high-ionization oxygen, neon, and sulfur indicate a source of far-UV photons 
in the hidden nucleus. However, we find that the other 55 ± 13% of narrow-line 
FR II radio galaxies are weak at 15 /mi, contrary to single-population unification 
schemes. Most of these galaxies are also weak at 30 /im. Mid-IR weak radio galax- 
ies may constitute a separate population of nonthermal, jet-dominated sources 
with low accretion power. 

Subject headings: galaxies: active, galaxies: quasars, galaxies: jets, infrared: 
galaxies 
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1. Unification of Quasars and Radio Galaxies 

The nature of the energy source in active galactic nuclei (AGNs) is a fundamental 
problem. The basic model attributes the large luminosity of these systems to gravitational 
energy release in an accretion disk around a supermassive black hole. A jet may be driven by 
magnetic fields threading the disk (Blandford & Payne 1982). The black hole spin energy 
may also be tapped and converted into electromagnetic Poynting flux and particles in a 
relativistic jet (Blandford & Znajek 1977; Punsly & Coroniti 1990; Meier 1999; De Villiers 
et al. 2004). 

Extragalactic radio sources are categorized by their morphology as either of two types 
(Fanaroff & Riley 1974). Fanaroff- Riley (FR) type I sources are edge-darkened, while FR 
lis are edge-brightened. The different morphology of FR Is indicates that they are not 
related to FR lis by orientation. FR Is also have lower radio luminosities than FR lis 
for a given host galaxy luminosity (Owen & Ledlow 1994; Bicknell 1996) and most have 
low-ionization nuclear emission region (LINER) spectra (Hine & Longair 1979). However, 
not all FR Is can be characterized by low accretion power (Whysong & Antonucci 2004; 
Cao & Rawlings 2004). The present paper focuses on FR lis, which contain powerful jets 
with bright terminal hot spots and lobes. Furthermore, we count broad-line radio galaxies 
(BLRGs) as low-luminosity quasars. 

Quasars and narrow-line radio galaxies (NLRGs) may be unified by orientation-dependent 
obscuration. Radio galaxies are thought to host quasar nuclei that are obscured by circum- 
nuclear dusty tori aligned with the radio jets (Antonucci 1984) . Unification of radio galaxies 
and quasars can therefore explain the lack of quasars viewed at large angles to the radio axis 
(Barthel 1989). The percentage of high-redshift radio galaxies (60% of the 3CRR FR II 
sample at z > 0.5) would then indicate a torus covering fraction of ~ 0.6. 

However, there appears to be a discrepancy between the redshift distributions of quasars 
and radio galaxies at z < 0.5, with a factor of ~ 4 more narrow-line radio galaxies than 
quasars (Singal 1993). Furthermore, the median projected linear size of these 'excess' radio 
galaxies is smaller than expected for quasars seen in the sky plane (Singal 1993; Whysong 
2005). The unification hypothesis may be modified to include a second population of lower 
luminosity, low-excitation FR II radio galaxies (Wall & Jackson 1997; Grimes, Rawlings, 
& Willott 2004). Alternatively, it has been argued that the torus covering fraction may 
increase with decreasing radio luminosity (Lawrence 1991). 

The unification hypothesis has been qualitatively confirmed by spectropolarimetry of 
radio galaxies, many of which have been shown to have highly polarized broad emission lines 
and blue continuum, scattered from material which has a direct view of the active galactic 
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nucleus (Cimatti et al. 1997; Cohen et al. 1999). Of particular note are the original 
discovery of highly polarized broad Ha from the hidden quasar nucleus in 3C 234 (Antonucci 
1984), and the discovery of highly polarized broad Ha in the spectrum of the powerful radio 
galaxy Cygnus A (Ogle et al. 1997). However, this method of detecting hidden quasars 
relies on an appropriately placed scattering region to view the otherwise hidden nucleus. 
Such a region is not guaranteed to exist for all radio galaxies, and thus spectropolarimetry 
can easily yield false negatives. Polarimetry is also ineffective at determining the luminosity 
of the hidden nucleus, since the scattering efficiency is usually unknown. 

Another way to search for hidden quasar nuclei is to observe radio galaxies in the mid-IR. 
If the unification hypothesis is correct, the dusty torus should serve as a crude calorimeter of 
the central engine (Meisenheimer et al. 2001; Siebenmorgen et al. 2004; Haas et al. 2004; 
Whysong & Antonucci 2004). Optical, UV, and X-ray photons from the quasar nucleus 
are absorbed by dust in the torus and the energy is re-emitted in the thermal infrared. 
This explains why blue, UV color-selected quasars emit 10-50% of their luminosity in the IR 
(Sanders et al 1989; Haas et al. 2000). There appears to be no connection between the bulk 
of this IR emission and nonthermal radio emission, except in core-dominated radio sources 
such as blazars. Observations of matched 3CR quasars and radio galaxies by ISO indicate 
similar IR luminosities, consistent with the unification picture (Meisenheimer et al. 2001; 
Haas et al. 2004). However, differences in 24 /iin/ 70 fxm color may indicate that mid-IR 
emission from the torus is anisotropic by a factor of < 3 (Shi et al. 2005). 

We present Spitzer observations of a sample of 42 FR II radio galaxies and quasars 
selected from the 3CRR survey. The goals are to search for mid-IR emission from hidden 
quasar nuclei and test the ubiquity of the unification hypothesis. The Spitzer Infrared 
Spectrograph (IRS) combines the advantages of unprecedented sensitivity from 5-36.5 /iin to 
measure the mid-IR continuum and spectral resolution to measure high ionization emission 
lines powered by hidden AGNs. In the current paper, we present evidence for hidden quasar 
nuclei based on mid-IR photometry extracted from the IRS spectra. We examine in detail the 
spectra of the subset of 14 mid-IR luminous radio galaxies which appear to contain hidden 
quasar nuclei. Spectra of the quasars and mid-IR weak radio galaxies and a statistical study 
of the complete sample will be presented in separate papers. 

2. Sample 

We begin by selecting a well-defined, radio flux-limited and redshift-limited sample of 
55 radio galaxies and quasars from the 3CRR catalog (Laing, Riley, & Longair 1983). We 
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include all 3CRR sources with FR II radio morphology, a flux of Sns > 16.4 Jy 1 at 178 MHz, 
and a redshift of z < 1. The original 3CRR catalog has a flux limit of 10 Jy at 178 MHz, 
is restricted to northern declinations (5 > 10°), and has galactic latitude \b\ > 10°. It is the 
canonical low-frequency selected catalog of bright radio sources, has optical identifications 
and redshifts for all entries, and has been extensively observed in most wavebands. 

We select only sources with FR II radio morphology. We verify or update the FR 
classification of all sources by inspection of the latest published radio maps. Compact, steep- 
spectrum sources (CSSs: 3C 48, 138, 147, 286, and 309.1) with radio major axis D < 10 kpc 
(Fanti et al. 1985) are excluded from the sample because they may constitute a class of 
young or frustrated radio sources. Here and throughout this paper, we assume a cosmology 
with H = 70 km s _1 Mpc -1 , Q m = 0.3, and Q\ = 0.7. Size and morphology indicate that 
CSSs are not related to FR lis by orientation. 

It is essential for our unification studies that we select a sample based on isotropic 
radio lobe flux, and not on optical or IR properties, so that it is unbiased by orientation- 
dependent selection effects. In particular, our sample includes no blazars. No sources make 
the flux limit only because of beamed emission from the core of the radio jet. Our sample 
includes quasars as well as radio galaxies, and we use the quasar subsample as a control. We 
aim to determine whether and which narrow-line FR II radio galaxies have mid-IR power 
comparable to quasars or broad-line radio galaxies of similar radio lobe flux and redshift. 

The 42/55 sources in our sample which we have observed with Spitzer, or which have 
Spitzer data in the public archive are listed in Tables 1 and 2. The 25 mid-IR luminous 
sources with uL u (15fim) > 8 x 10 43 erg s _1 (14 NLRGs and 11 quasars or BLRGs) are listed 
in Table 1, and the 17 mid-IR weak galaxies with uL u (15fim) < 8 x 10 43 erg s _1 are listed 
in Table 2. The reason for this particular division is explained below. 

Optical source classifications are based on emission line properties. Type 1 sources have 
directly visible broad emission lines (quasars and BLRGs), and type 2 sources (NLRGs) 
do not. The NLRGs are further classified using their forbidden emission lines (Jackson & 
Rawlings 1997; Wilott et al. 1999) 2 . High-excitation galaxies (HEGs) are defined to have [O 
ill] A5007 equivalent widths of > 10 A and [O ill] A5007/[O n] A3727 > 1. The sources which 
do not meet these criteria are classified as low-excitation galaxies (LEGs). The equivalent 
width criterion ensures that [O ill] is measurable in moderate S/N spectra. However, it 



1 Thc radio flux limit is 15 Jy using Laing, Riley, & Longair (1983) flux values and 16.4 Jy on the standard 
Baars (1977) scale. 

2 Updated optical classifications are available at http://www-astro.physics.ox.ac.uk/~cjw/3crr/3crr.html. 
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remains to be seen whether some sources with low [O ill] equivalent width might have [O 
ill] A5007/[O n] A3727 > 1. In addition, we caution that [O n] and [O ill] may be subject 
to differing amounts of extinction. 

3. Observations 

We observed the sources in our sample with the Infrared Spectrograph (IRS) on the 
Spitzer Space Telescope (Houck et al. 2004; Werner et al. 2004). We used the low-resolution 
(R ~ 64 — 128) modules Short-Low (SL) and Long-Low (LL) for accurate spectrophotometry 
over the wavelength range 5-36.5 /xm. Wavelengths 36.5-40 fim are unusable because of 
low-S/N and 2nd-order bleed-through caused by filter delamination in LL 1st order (LL1). 
The absolute and relative flux accuracies of IRS are generally better than 10% and 4%, 
respectively, as judged from observations of bright standard stars. However, additional low- 
level instrumental artifacts may become important for faint sources. 

We used IRS in standard flux-staring mode, for 2 cycles at 2 nod positions in each of the 
modules SL1, SL2 (SL 1st and 2nd order), and LL2 (LL 2nd order). We executed 1 cycle at 
2 nod positions for LL1, which covers the 20-36.5 /im range. A typical observation includes 
240 s of on-source exposure time in each of SL1, SL2, and LL1, and 480 s in LL2, for a total 
of 2000 s per target (including overhead). Archived IRS data are used for 14 sources which 
were observed partly or in full by other investigators (with similar or longer exposure times). 

Nod or off-slit observations were subtracted to remove foreground emission from the tele- 
scope, zodiacal light, and interstellar medium. Spectra were then extracted from the Basic- 
Calibrated Datasets (BCDs), using the Spitzer IRS Custom Extraction (SPICE 3 version 1.1) 
software and standard tapered extraction windows. The extraction window full-widths are 
proportional to wavelength in each order to match the diffraction-limited telescope point- 
spread function (SL2: 7"2 at 6 /im, SL1: 14'/4 at 12 /im, LL2: 21'/7 at 16 fim, LL1: 36"6 at 
27 pm). We rebinned portions of the spectra of 3C 55, 172, 220.1, 244.1, 263.1, 280, and 
330 by factors of 4-8 in order to improve the S/N at short wavelengths. Spectral orders were 
trimmed at the edges and merged to produce final spectra. 

The SL and LL slits have widths of 3'.'7 and 10'.'6, respectively. Standard point-source flux 
calibrations (version 12.0) were applied to correct for slit and aperture losses and convert 
the spectra from electron s _1 to Jy. In most cases, fluxes match to < 15% across order 
boundaries, consistent with a point source that is well-centered in all of the slits. However, 
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in 5 cases (3C 192, 216, 220.1, 380, and 381) SL2 fluxes are larger by 17-35% relative to the 
other orders. Assuming that these mismatches owe to variable slit-loss caused by pointing 
errors, the orders with low flux are adjusted upward to match the orders with high flux. 
Order mismatches may alternatively be an indication of extended mid-IR emission. 

The results for a few sources with nearby neighbors in the slit should be viewed with 
caution. In the case of 3C 310, a nearby companion galaxy (to the east) may contribute 
a significant fraction of the flux (< 50%) in the LL1 slit. Similarly, a nearby source may 
potentially contribute to the LL spectra of 3C 438 (which is, however undetected at 15 /im). 
The SL2 spectrum of 3C 388 may be weakly affected by flux from a nearby star on the slit 
(< 20%). The northern component of the double nucleus in 3C 401 falls outside of the SL 
slits, but falls inside the LL slit used to measure the 15 /im flux. 

3.1. Mid-Infrared and Radio Luminosities 

We measure the mean 6.5 — 7.5 /im and 13.0 — 17.0 /im flux densities F v (7 and 15 /im, 
rest) of each target (Tables 1 & 2). All SpitzerHux densities in this paper are in observed units 
at a constant rest-frame wavelength defined by A res t = A bs/(1 + z), where z is the redshift 
measured from optical emission lines and cataloged in the NASA Extragalactic Database 
(NED 4 ). This avoids any complication from potentially large cosmological K-corrections that 
could otherwise be introduced by a steep IR continuum slope or redshifted silicate absorption 
features. 

We choose to measure the mid-IR flux at 7 and 15 /tm to avoid the 9.7 /im trough and 
the deepest part of the 18 /im silicate absorption trough. We exclude the 14.0-14.5 /im and 
15.3-15.8 /im wavelength regions from our photometry, to avoid emission from Ne V and Ne 
ill. The 7 and 15 /im bands are within the Spitzer IRS bandpass for redshifts z < 1.28. 
However, the archival LL data for two quasars (3C 254 and 275.1) are not yet public. We 
extrapolate their SL spectra to obtain F v (15 /im, rest) using F u (7 /im, rest) and the observed 
(relatively line-free) 5-7 /im spectral index. 

Radio luminosities uL u (178 MHz, rest) are estimated from the observed 178 MHz fluxes 
and K-corrected using the 178-750 MHz radio spectral index (Laing, Riley, & Longair 1983). 
The sources in our sample display a large range of nearly 3 orders of magnitude in mid-IR to 
radio luminosity: vL u (15 /im)/VL„(178 MHz) = 0.8 — 680 (Fig.l). This quantity is thought 
to reflect the relative importance of accretion luminosity and jet kinetic power dissipation. 
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However, different size and time scales are probed by the radio (10 kpc-1 Mpc) and mid-IR 
(0.1-100 pc), and the radio power may be sensitive to differences in environmental conditions. 

For the purpose of studying quasar and radio galaxy unification, it is natural to divide 
the sample into the mid-IR luminous NLRGs which emit as powerfully as quasars or BLRGs, 
and the mid-IR weak NLRGs that do not. We adopt an empirical dividing line of uL u (15 
jum) > 8 x 10 43 erg s" 1 to separate hidden quasars from mid-IR weak radio galaxies. The 
cutoff is set at 1/2 the luminosity of the mid-IR weakest BLRG (3C 219) to allow for some 
degree of anisotropy at 15 pm. Fourteen NLRGs satisfy our criterion and are thus likely to 
contain hidden quasar or BLRG nuclei (Table 1). Notably, all of these NLRGs are optically 
classified as HEGs. 

The 17 mid-IR weak NLRGs with i/L„(15 /mi) < 8 x 10 43 erg s' 1 (Table 2) have 
mixed optical classifications, including both HEGs and LEGs. These sources have lower 
S/N mid-IR spectra, which will be considered in detail in a later paper. Six mid-IR weak 
NLRGs (including 2 HEGs and 4 LEGs) are undetected by Spitzer at 15 /im, and one is also 
undetected at 7 /im. 

3.2. Hidden Quasar Spectra 

3.2.1. Continuum Emission 

We now present Spitzer spectra of the 14 mid-IR luminous NLRGs that ostensibly con- 
tain hidden quasar nuclei (Figs. 2-4). We also plot the spectral energy distributions (SEDs) 
of the sources with published near-IR photometry (Fig. 5). The collected photometric data 
were measured in the J, H, K, L', and M wavelength bands from the ground (Lilly & Longair 
1984; Lilly, Longair, & Miller 1985; Simpson, Rawlings, & Lacy 1999; Simpson & Rawlings 
2000). The photometric apertures range in size from 3 — 11", with preference given to the 
apertures that most closely match the Spitzer SL slit width. Where available, the ground- 
based L' and M-band photometry agrees with Spitzer spectrophotometry remarkably well. 
There is no indication of variability over the time span of 20 yr. 

Four of the low-redshift NLRGs (3C 33, 234, 381, and 452) have broad peaks in their vL v 
spectra (and SEDs) at 1.5 — 2.5 x 10 13 Hz (12-20 /im). A maximum and spectral curvature 
near 20 /xm are also suggestive of broad peaks in the Spitzer spectra of 3C 55, 244.1, 265, 
and 330. The large amplitude (~ 0.5 — 1.0 dex) of the mid-IR bump (Fig. 5) excludes a 
large contribution of synchrotron emission to the mid-IR continuum of most sources. This 
is not surprising if the equatorial plane of the dusty torus is roughly perpendicular to the 
radio jet, such that jet emission is beamed away. The high redshifts of the NLRGs 3C 172, 
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220.1, 263.1, 268.1, and 280 preclude the identification of a mid-IR bump in the SEDs of 
these sources. The unusually flat, blue SED of 3C 433 may indicate a quasar viewed at low 
inclination (Section 3.2.2). 

We attribute the mid-IR continuum bump visible in most sources to thermal emission 
from warm or hot dust. Fitting the mid-IR peak with a single-temperature blackbody model 
indicates dust with a temperature of 210 — 225 ± 0.5 K (Fig. 5). While this temperature 
characterizes the peak of the mid-IR SED, hotter dust must also be present. At frequencies 
greater than the peak of the SED (2.0 — 7.5 x 10 13 Hz), the continuum emission of most 
sources can be characterized using a power law with spectral index a = 1.1 — 2.1 (Table 1 & 
Fig. 6). This emission likely comes from a continuous distribution of dust temperature. We 
measure the spectral index between 7 and 15 /im, avoiding the 9.7 /im and 18 /im silicate 
absorption troughs. The most blue and apparently hottest mid-IR luminous NLRG is 3C 
265, while the most red and coolest are 3C 55 and 3C 268.1 (Fig. 6). In comparison, some 
mid-IR weak sources such as 3C 310 and 3C 388 are quite blue (a ~ —0.1 — h0.7), indicating 
a large contribution of starlight from the host galaxy to the 7 /im continuum. 

The near-IR continuum shifts into the Spitzer IRS passband for the highest redshift 
(z > 0.7) sources. The spectra of the NLRGs 3C 55 and 3C 265 steepen above 7.5 x 10 13 
Hz (below 4 pm). Fitting these spectral turnovers with single-temperature blackbodies, we 
find emission from hot dust with temperatures of 520 ± 10 K and 660 ± 10 K, respectively. 
Altogether, the mid-IR luminous radio galaxies in our sample show emission from dust 
with temperatures distributed in the range 210-660 K. Hotter temperature dust (up to the 
sublimation temperature) may be present but not visible for radio galaxy tori viewed at high 
inclination (Pier & Krolik 1992). 

Extinction by cold foreground dust in the host galaxy may also affect the spectral 
index. For Galactic-type dust, A(7, 15, 35 /im)/A(V) = (0.020,0.015,0.004) (Mathis 2000). 
An extinction of A(V) = 100 would steepen the 7-15 /im spectral index by 5a = 0.6 (Fig. 
6a). The observed range in spectral index for the mid-IR luminous NLRGs is 5a = 1.0, 
corresponding to A(V, 7, 15, 35 /im) = (167, 3.3, 2.5, 0.6) mag. Thus if reddening by a cold 
foreground dust screen accounted entirely for the range in mid-IR slope, the mid-IR emission 
could be anisotropic by factors of /a(7, 15, 35 /im) ~ (22, 10, 1.3). However, these are upper 
limits since variations in the spectral index are also controlled by the physical temperature 
distribution of the visible dust. 

The SEDs of several sources (3C 33, 55, 172, 265, and 452) have upturns at short 
wavelengths, which we attribute to stellar emission from the host galaxy (Fig. 5). The 
wavelength of the upturn (1-5 /im) is an indicator of the relative strength of the mid-IR bump 
seen from our direction vs. host galaxy light, occurring at shorter wavelength for sources 
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with a stronger mid-IR bump. This may have important consequences for understanding 
the K-z Hubble diagram for 3C radio galaxies, for which it has been argued that AGNs 
contribute a negligible fraction of the K-band flux (e.g., Simpson & Rawlings 2000). This 
may be incorrect for a few of the most luminous mid-IR sources in our sample, including 3C 
234 and 3C 280 where there appears to be much emission from hot dust in the K band. 

Detailed spectral modeling, combined with radio orientation indicators, promises to 
further characterize the temperature distribution, optical depth, and inclination of the dusty 
torus that is thought to produce most of the mid-IR emission from hidden quasar nuclei. 
Such an analysis is, however, outside the scope of the present paper. 

3.2.2. Silicate Absorption 

The silicate absorption trough at 9.7 /im is detected in 12/14 of the mid-IR luminous 
NLRG spectra (Table 3 & Figs. 2-4). The equivalent width EW9.7 and apparent optical depth 
T9.7 are measured relative to a local continuum fit to either side of the trough, indicated in 
Figures 2-4. The optical depth is averaged over the trough bottom (rest 9.2-10.2 /jm) to 
improve the S/N. It should be kept in mind that the apparent r 9 . 7 is just a convenient 
parameterization of (and lower limit to) the total optical depth since there must also be 
broad-band silicate absorption of the adjacent continuum. 

The apparent silicate optical depths are small (r 9 . 7 = 0.02—0.2), for all but 3C 55 and 3C 
433. If attributed to foreground dust screens, this would indicate optical extinction of only 
Ay = 0.2 — 5.1 mag (Fig. 6b), assuming a Galactic extinction law with Ay/r 9 . 7 = 12.3 — 25.6 
mag (Rieke & Lebofsky 1985; Draine & Lee 1984). The extinction values are clearly 
underestimated since they imply that the hidden nuclei in 3C 234, 265, 381, and 452, which 
have T9.7 < 0.1, should be reddened but directly visible at Ha. The same discrepancy between 
T9.7 and estimates of extinction at shorter wavelengths is seen for the hidden quasar nucleus 
in Cygnus A, and attributed to a radial gradient in torus dust temperature (Imanishi & 
Ueno 2000). The observed range of T9.7 may correspond to a range of equatorial silicate 
dust column densities in the torus, or alternatively a range of viewing angles. In this regard, 
more detailed modeling of the torus, including its geometric and temperature structure is 
clearly called for. 

Filling-in of the silicate troughs by silicate emission from the torus or narrow-line region 
(NLR) may also reduce the apparent silicate optical depths in some sources. This is predicted 
for an optically thick torus viewed at an intermediate or face-on inclination (Pier & Krolik 
1992). Recently, strong silicate emission features were detected by Spitzer in several radio- 
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loud (3C) and radio quiet (PG) quasars (Siebenmorgen et al. 2005; Hao et al. 2005). 
The failure of previous attempts to observe this feature inspired torus models with large 
dust grain size (Laor & Draine 1993) or a spatial distribution of optically thick clumps 
(Nenkova, Ivezic, & Elitzur 2002). However, it appears that past non-detections owe to 
inadequate wavelength coverage to determine the underlying continuum. 

The NLRGs 3C 55 and 3C 433 have significantly deeper silicate troughs than other 
NLRGs, with T9.7 = 0.9 and 0.7, respectively (Fig. 6b). We suggest that their active nuclei 
and tori are absorbed by an additional (kpc-scale) cold dust screen in the host galaxy. As 
noted above, the NLRG 3C 433 is unusual in having a flat, blue continuum (similar to some 
of the quasars in our sample). A blue mid-IR spectrum is not necessarily at odds with deep 
silicate absorption features. It can be understood if this is a quasar viewed at low inclination 
to the jet and torus axes, but through an (Ay ~ 10) cold dust screen. This amount of 
extinction would result in very little reddening at 7-15 /im (5a = 0.05 — 0.14), but would be 
sufficient to create the deep 9.7 /im trough (Fig. 6b) and would obscure any optical broad 
lines. 

The NLRG 3C 433 is also unique in having the only unambiguously detected 18 /iin 
silicate trough, with equivalent width EWi 8 = 0.42 ± 0.01 fim and apparent optical depth 
Ti8 = 0.07 ± 0.03 (averaged over 17-19 /xm). The ratio of rig to T9.7 apparent silicate trough 
depths is 0.10 ± 0.04, consistent with a 0.11 ratio for Galactic-type silicate dust (Draine & 
Lee 1984). We do not see the full 18 /zm silicate trough in the spectrum of 3C 55 because 
of inadequate rest- wavelength coverage. 

3.2.3. Forbidden Emission Lines 

All of the mid-IR luminous NLRGs with high S/N spectra have forbidden emission lines 
from highly ionized metals, including [O iv] A25.89 /im, [Ne 11] A12.81, [Ne ill] A15.55, [Ne 
v] A14.3, [Ne v] A24.31, [Ne vi] A7.65, [S in] A18.71, [S in] A33.48, and [S iv] A10.51 (Figs. 
2-4). We measure the line flux and rest equivalent width of each emission line relative to 
the local continuum level (Table 4). Formal uncertainties are computed from the noise in 
the continuum to either side of the line. Upper limits are estimated for undetected emission 
lines, assuming they are unresolved. 

The large range of ionization states (especially high-ionization Ne V, Ne VI, and S iv) 
indicates photoionization by a hidden source of far-UV photons (Alexander et al. 1999; 
Sturm et al. 2002; Armus et al. 2004), e.g. a quasar nucleus. Low critical densities in the 
range 10 3 — 10 6 cm -3 (Alexander et al. 1999) indicate that the forbidden lines arise in the 
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NLR. There could plausibly be contributions from starburst emission to the lower-ionization 
emission lines such as [Ne n]. For 3C 55 and 3C 433, the [S iv] A10.51 line has a relatively 
large flux even though it sits at the bottom of a deep silicate trough. This line must then 
arise from a region not heavily obscured by dust, such as the extended NLR. The resolving 
power of IRS is insufficient to measure the intrinsic emission line widths, which are therefore 
< 4700 km s _1 . 

In order to assess the ionization level of the emission line regions, we compute several 
emission line ratios (Fig. 7). In particular, the [O lv]/[Ne n] and [Ne v]/[Ne n] ratios can be 
used as diagnostics of the relative contributions of AGN and starburst emission to the mid- 
IR emission line spectra of galaxies (Genzel et al. 1998; Sturm et al. 2002). In the sources 
where O IV , Ne V, and Ne II emission are all detected (3C 33, 234, 381, and 433), we find [O 
iv] A25.89 jum / [Ne il] A12.81 /mi > 1.0 (Fig. 7a) and [Ne v] A14.3/ [Ne il] A12.81 /im > 0.5 
(Fig. 7c), indicating a > 50% AGN contribution to the emission line spectrum. The [Ne il] 
A12.81 /im line is relatively weak or undetected in the z > 0.2 sources, making it difficult to 
apply these diagnostics. However, the large EWs of the [Ne vi] or [S iv] lines in 3C 55, 265, 
and 330 indicate that the emission line spectra of these sources are also AGN-dominated. 



3.2.4- Molecular Emission 

Polycyclic aromatic hydrocarbon (PAH) emission features are commonly seen in star- 
forming regions, starbursts, and starburst-dominated ULIRGs (e.g., Genzel et al. 1998; 
Armus et al. 2004). The only PAH feature we detect is the weak 11.3 /im line in the 
spectrum of 3C 33, with a flux of 2.3 ± 0.3 x 10~ 14 erg s _1 cm~ 2 and equivalent width of 
0.009 ± 0.003 /im (Fig. 2). We do not detect a 6.2 /mi PAH feature in 3C 33 (EW< 0.8 
/mi) or any of the other mid-IR luminous NLRGs, though this spectral region is generally 
noisier. If present, we could not cleanly resolve the 7.7 /im PAH feature from the adjacent 
[Ne Vi] line, nor the 12.7 /mi PAH feature from the adjacent [Ne il] line. Regardless, we do 
not see any hint of these PAH features in any source, suggesting EW<< 0.1 /im. Therefore 
in most cases, neither the 7.7 or 11.3 /im PAH features can have a significant impact on the 
measurement of the silicate trough. 

The general lack of PAH features is a strong indication that the primary power source 
in mid-IR luminous radio galaxies is accretion power, not hot stars. It is likely that PAHs 
are destroyed in the torus, which is exposed to intense X-ray emission from the AGN (Voit 
1991). The weak PAH emission that is present in 3C 33 may arise in star-forming regions 
shielded from the AGN. 
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We detect the H 2 0-0 S(3) 9.67 fim and H 2 0-0 S(l) 17.03 jum pure rotational emission 
lines of molecular hydrogen (at the 3a level) only in the spectrum of the NLRG 3C 433 (Fig. 
2). The line fluxes are 0.6 ±0.2 x 10~ 14 and 1.4 ±0.4 x 10~ 14 erg s _1 cm -2 respectively 
(EW = 0.008 and 0.013 /im). The location of the 9.67 fim line at the bottom the deep 
9.7 fim silicate trough may indicate that the H2 emission region is exterior to the obscuring 
dust screen. The H 2 emission lines can be produced in warm molecular gas heated either by 
shocks or X-ray photons from the AGN (Rigopoulou et al. 2002). The ratio of S(3)/S(l) 
line fluxes is 0.5 ± 0.2, which indicates warm H 2 with an excitation temperature of 300 ± 30 
K. We estimate a warm H 2 mass of roughly 2 x 1O 8 M , assuming a Boltzmann distribution 
of rotational level populations, an unresolved source, and negligible mid-IR extinction. 



4. Discussion 

4.1. Hidden Quasar Nuclei 

The high mid-IR luminosities vL u (15 /im) = 10 44 - 10 46 erg s" 1 of 45 ± 12% (14/31) of 
the FR II NLRGs in our sample are consistent with hidden quasar or BLRG nuclei. In fact, 
such copious hot dust emission directly requires a hidden source of quasar-like luminosity to 
power it. Including the 11 quasars and BLRGs seen directly, we find that at least 60 ± 12% 
(25/42) of 3CRR FR II sources at z < 1 with S l78 > 16.4 Jy contain powerful AGNs. 

The percentage of mid-IR luminous AGNs obscured by dust and therefore classified as 
NLRGs is 56 ± 15% (14/25), corresponding to a mean torus covering fraction of 0.56 and 
mean torus opening half-angle of 55 ±11°. (If the 8 mid-IR weak HEGs are counted as highly 
obscured quasars, then the mean torus covering fraction increases to 0.67 ± 0.14.) Both of 
these numbers are consistent with the estimated 60% mean torus covering fraction required 
to unify z = 0.5 — 1.0 radio galaxies and quasars (Barthel 1989). 

The receding torus model (Lawrence 1991) predicts a larger torus covering fraction 
for low luminosity sources. If so, we would expect a larger mean mid-IR/radio ratio and a 
smaller type 1 fraction for low-redshift than for high-redshift sources. The fraction of type 1 
mid-IR luminous sources is 3/9 (0.3 ±0.2) at z < 0.5 vs. 8/16 (0.5 ±0.2) at z > 0.5. Clearly, 
a larger sample is needed to put meaningful constraints on the variation of torus covering 
fraction with redshift or luminosity. 

For the galaxies that contain a powerful mid-IR source, we find that there are other 
indications of a hidden AGN. The high- ionization, mid-IR forbidden lines such as [Ne v], 
even more-so than the strong optical [O ill] lines, are telltale signatures of a non-stellar source 
of FUV photons in the radio galaxies that show them. Because of the lower extinction in 
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the mid-IR relative to the optical, it is likely that we can see these lines closer to the nucleus 
than optical narrow lines such as [0 ill] (Haas et al. 2005). 

At least three of the mid-IR luminous NLRGs in our sample are known to have highly 
polarized broad emission lines. The NLRG 3C 234 has quasar-like mid-IR luminosity and a 
highly polarized broad Ha line (Antonucci 1984). The NLRG 3C 265 has quasar luminosity, 
highly polarized UV flux, and a highly polarized broad Mg II line (Tran et al. 1998). The 
low-redshift NLRG 3C 33 has mid-IR luminosity comparable to the BLRG 3C 219, and 
highly polarized broad Ha and B./3 (Cohen et al. 1999). We have an ongoing program to 
obtain optical spectropolarimetry of the rest of the sources in our sample. While the radio 
galaxies with highly polarized broad lines are known to contain hidden type-1 AGNs, there 
were previously no reliable measurements of the hidden AGN luminosities. Our mid-IR flux 
measurements yield rough calorimetric estimates of the hidden AGN luminosities, subject 
to uncertainties in SED, torus covering fraction, and any mid-IR anisotropy. 



4.2. Mid-IR Weak Radio Galaxies 

The majority of radio galaxies in our sample (17/31 or 55 ± 13%) are relatively weak 
mid-IR sources. It is possible that some of the AGNs are highly obscured even at 15 /zm 
because they are viewed through a very large dust column. A mid-IR source obscured by a 
nearly Compton-thick (r c = 0.7) disk with Galactic dust/gas ratio would have an equatorial 
extinction of ~ 5 mag (factor of 100) at 15 fim. Even in this case, any mid-IR emission 
above the disk might not be obscured. For example, there may be a contribution from dust 
in the narrow-line region (NLR), above the hole in the torus (e.g. NGC 1068, Galliano et al. 
2005; Bock et al. 2000), tending to make the mid-IR emission more isotropic. 

For galaxies with redshift z < 0.22, Spitzer can measure the flux at A = 30 jum (rest), 
which should be more isotropic and less subject to extinction than the 15 /im flux. Nev- 
ertheless, 9/11 mid-IR weak galaxies in this redshift range are also weak at 30 /im, with 
vL v (30 fim) < 8 x 10 43 erg s _1 . The two exceptions are 3C 61.1 and 3C 123, with uL u (30 
/im) = 1.27 ± 0.08 and 1.4 ± 0.2 x 10 44 erg s -1 , respectively. In comparison 3C 452, the 
weakest mid-IR luminous NLRG in this redshift range, has uL„(30 /im) = 8.63 ±0.09 x 10 43 
erg s -1 . Therefore, reclassifying the NLRGs by their luminosity at 30 //m would only gain us 
an additional 2 mid-IR luminous sources. This leads us to believe that most of the mid-IR 
weak radio galaxies truly lack a powerful accretion disk. Relatively low accretion power 
suggests, but does not prove, that some FR II jets may be driven by radiatively inefficient 
accretion flows or black hole spin-energy (Begelman, Blandford, & Rees 1984; Meier 1999). 
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As we mentioned above, roughly half (9/17) of the mid-IR weak NLRGs are are classified 
as LEGs with weak optical [0 ill] emission. The [O ill] emission in these sources may be 
weak because there is no strong source of UV photons to power the NLR. Qualitatively 
similar conclusions have been drawn by other investigators (e.g., Hine & Longair 1979; 
Chiaberge et al. 2000; Grimes, Rawlings, & Willott 2004). Alternatively, the NLR may be 
partly or completely obscured in LEGs (Haas et al. 2005). It will be important to make 
a more quantitative assessment of the optical and IR emission line strengths, to evaluate 
the extinction and determine what UV luminosity is necessary to power the NLR in mid-IR 
weak NLRGs. 



4.3. Radio Properties and Unification 

One of the major motivations for the radio galaxy and quasar unification theory is 
the deficit of lobe- dominant vs. core-dominant FR II quasars (Barthel 1989). Relativistic 
beaming models predict that there should be relatively more sources where the radio jet is 
beamed away and the high frequency radio core is weak. Core fluxes at 5 GHz are tabulated 
for the 3CRR catalog by Laing et al. (2003) 5 . We identify the mid-IR luminous NLRGs in 
our sample with the missing lobe-dominant quasar population (Fig. 8). Their median core 
to lobe ratio is R c = z/L^core, 5 GHz)/z/L^(178 MHz) = 5, while the median R c = 180 for 
the quasars. Conversely, all of the mid-IR luminous sources with R c > 100 are classified as 
quasars. We will perform a more detailed statistical analysis when the Spitzer observations 
of our full sample are complete. 

For most mid-IR weak NLRGs, we reject the possibility that the AGN and radio core 
are viewed in an 'off' state while the radio lobes are still active. A 5 GHz core is detected 
in 13/17 of the mid-IR weak radio galaxies (Fig. 8) and 13/14 of the mid-IR luminous 
galaxies. Furthermore, mid-IR weak and mid-IR luminous NLRGs have comparable core to 
total luminosity ratios of R c — 1 — 100. Deeper 5 GHz observations of the 5 non-detected 
cores (in 3C 28, 153, 172, 315, and 319) will be necessary to determine whether or not they 
are in an off state (e.g., R c < 0.1). 

We find that FR II radio morphology is not a reliable predictor of nuclear mid-IR 
luminosity for radio galaxies. Contrary to the simple unification paradigm, not all narrow- 
line FR II galaxies host nuclei as powerful as quasars with matched radio lobe luminosity. 
Unification theories must be modified to account for an additional population of mid-IR 
weak radio galaxies. 



5 The online 3CRR catalog is available at http://www.3crr.dyndns.org. 
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Both intrinsic jet power and interaction with the interstellar and intergalactic medium 
are likely to be important for determining radio morphology. The break luminosity between 
FR I and FR II radio sources is found to increase with host galaxy optical luminosity (Owen & 
Ledlow 1994; Bicknell 1996). The existence of radio sources with hybrid FR I/II morphology 
also points to the importance of environmental effects in determining radio morphology 
(Gopal-Krishna & Wiita 2000; Gawronski et al. 2005). Furthermore, most FR I radio jets 
are one-sided, relativistic, and narrowly collimated on sub-parsec scales, just like FR lis, and 
decollimate only on kpc scales (e.g. M87, Junor, Biretta, & Livio 1999). 

Contrary to a common misconception, not all FR I sources have radiatively inefficient 
nuclei. For example, Centaurus A is persuasively argued to have a powerful hidden AGN 
(Whysong & Antonucci 2004), and the BLRG 3C 120 is a well-known FR I source. Deep 
VLA observations of the optically luminous, 'radio-quiet' quasar E1821+643 demonstrate 
that it has an FR I radio morphology (Blundell & Rawlings 2001). These and other cases 
(Antonucci 2001) demonstrate that many powerful AGNs are FR Is. 

The observed variation in radio morphology and a wide range in AGN radio-loudness 
(Kellerman et al. 1989) do not necessarily require a weak coupling between jet power 
and accretion power, but may demonstrate that multiple factors are at work. At least five 
parameters may be necessary to theoretically unify all AGN types: black hole mass, black 
hole spin, accretion rate, radiative efficiency, and viewing angle. There is still much work 
ahead before we completely understand how basic physical parameters regulate the activity 
of supermassive black hole systems. Models that tie jet production to accretion onto a 
spinning black hole are particularly promising (e.g., Meier 1999; Hawley & Krolik 2006). 

Much progress has been made in understanding the aspect-dependent appearance of 
AGN disks and jets, as a consequence of relativistic beaming and obscuration (Urry & 
Padovani 1995). Our Spitzer observations confirm that many FR II radio galaxies would 
appear as powerful quasars if viewed from an unobscured direction (e.g. along the radio 
axis). However, just as many FR II radio galaxies would not. Our Spitzer observations 
also demonstrate that powerful radio jets may be produced even by mid-IR weak AGN. A 
powerful, luminous accretion disk is not always necessary to produce a highly collimated, 
relativistic jet. 



5. Conclusions 



(1.) We report on a large Spitzer spectroscopic survey of 3CRR FR II radio sources 
with S 178 > 16.4 Jy at z < 1. We find strong mid-IR emission from 45 ± 12% (14/31) 
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of NLRGs, which have luminosities comparable to matched BLRGs and quasars. Other 
indicators including high-ionization mid-IR lines and highly polarized broad emission lines 
confirm that some of these sources contain hidden quasar or BLRG nuclei. This demonstrates 
the power of Spitzer IRS for unveiling hidden quasars and estimating their luminosities. 

(2.) We present Spitzer spectra of the 14 mid-IR luminous radio galaxies. In most 
cases, the mid-IR continuum bump from 3-30 /im can be produced by a distribution of hot 
dust with temperatures in the range 210-660 K. These high temperatures are most likely 
maintained by hidden AGNs. The silicate absorption trough at 9.7 fim has an apparent 
optical depth of r = — 0.2 in most cases, consistent with dust temperatures decreasing 
outward from the center of a dusty torus. Two sources, 3C 55 and 433, have deeper silicate 
troughs which may be produced by additional cool dust in the host galaxy. 

(3.) However, not all FR II radio galaxies emit strongly in the mid-IR. Contrary to 
single-population unification schemes, the majority of narrow-line radio galaxies in our sam- 
ple (17/31 or 55 ± 13%) have weak or undetected mid-IR emission compared to matched 
quasars and BLRGs, with uL u (15 /im) < 8 x 10 43 erg s -1 . For a few sources, this may possi- 
bly be the result of anisotropic torus emission viewed through a large column density of dust. 
However, it is likely that most of the weakest sources do not contain a powerful accretion 
disk. These may be truly nonthermal, jet-dominated AGNs, where the jet is powered by a 
radiatively inefficient accretion flow or black hole spin-energy rather than energy extracted 
from an accretion disk. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. We have also made use of the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with NASA. Support for this research was provided by NASA through an award 
issued by JPL/Caltech. We thank Dave Meier, Lee Armus, Bill Reach, and the anonymous 
referee for their helpful input and comments on the manuscript. 
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Table 1. Mid-IR Luminous Sources 
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a Optical spectral type (Lacy et al. 1999; Jackson & Rawlings 1997). 

b,c Flux densities (mJy) and 3cr upper limits at 7 fim (rest) and 15 /im (rest) from 
Spitzer IRS. 
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d Logarithm of luminosity (erg s x ) at 15 /im (rest). Values for 3C 254 and 3C 275.1 
are extrapolated from 7 /iin because the LL data are unavailable. 

c Spectral power law index for F u ~ u~ a from 7-15 /im (rest). 

f Radio flux density (Jy) at 178 MHz (observed) (Laing, Riley, & Longair 1983), 
multiplied by a factor of 1.09 to convert to the Baars (1977) standard flux scale. 

g The broad Ha line visible in total flux is entirely scattered (Antonucci 1984). 

h Spectral power law index for 3C 268.1 measured from 8-15 //m (rest). 
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Table 2. Mid-IR Weak Sources 
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0.7 ± 0.2 


43.40 


1.3 ± 0.4 


19.5 


315 


HEG 


0.1083 


0.97 ± 0.08 


1.9 ± 0.2 


43.01 


0.9 ± 0.2 


19.4 


388 


HEG 


0.0917 


0.98 ± 0.06 


0.84 ± 0.09 


42.66 


0.4 ± 0.2 


26.8 


436 


HEG 


0.2145 


0.65 ± 0.05 


1.5 ± 0.2 


43.52 


1.1 ± 0.2 


19.4 


438 


HEG 


0.290 


0.16 ± 0.04 


<0.45 


<43.27 




48.7 


28 


LEG 


0.1953 


0.45 ± 0.06 


<0.30 


<42.74 




17.8 


123 


LEG 


0.2177 


1.07 ± 0.08 


2.8 ± 0.4 


43.81 


0.7 ± 0.2 


206.0 


153 


LEG 


0.2769 


0.29 ± 0.05 


1.0 ± 0.2 


43.59 


1.7 ± 0.3 


16.7 


173.1 


LEG 


0.2921 


0.38 ± 0.04 


0.6 ± 0.1 


43.40 


0.6 ± 0.3 


16.8 


288 


LEG 


0.2460 


0.40 ± 0.05 


<0.60 


<43.25 




20.6 


310 


LEG 


0.0538 


0.81 ± 0.07 


0.73± 0.1 


41.98 


-0.1 ± 0.2 


60.1 


319 


LEG 


0.1920 


<0.15 


<0.27 


<42.68 




16.7 


326 


LEG 


0.0895 


0.67 ± 0.09 


<0.39 


<42.16 




22.2 


401 


LEG 


0.2011 


0.25 ± 0.05 


0.8± 0.2 


43.17 


1.5 ± 0.4 


22.8 



a -^See Table 1. 
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Table 3. Silicate Trough 



3C 


EW 9 . 7 a 


r 9.7 


b 


33 


-0.222 ± 0.007 


0.14 ± 


0.02 


55 


-1.51 ± 0.03 


0.9 ± 


0.1 


172 


> -0.4 






220.1 


> -0.3 






234 


-0.028 ± 0.006 


0.019 ± 


0.008 


244.1 


-0.23 ± 0.03 


0.18 ± 


0.08 


263.1 


-1.3 ± 0.2 






265 


-0.11 ± 0.03 


0.04 ± 


0.06 


268.1 


-0.4 ± 0.1 


0.2 ± 


0.5: 


280 


-0.23 ± 0.02 


0.16 ± 


0.09 


330 


-0.41 ± 0.06 


0.1 ± 


0.1 


381 


-0.19 ± 0.02 


0.07 ± 


0.02 


433 


-1.30 ± 0.01 


0.71 ± 


0.07 


452 


-0.196 ± 0.009 


0.10 ± 


0.02 


a The 


9.7 fim silicate 


trough 


(rest) 



equivalent width in /zm. 

b Apparent 9.7 /im silicate optical 
depth, averaged over 9.2-10.2 /im (rest). 
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Table 4. Emission Lines' 



op 
Ob 


[INe vij 


[S IVj 


["AT.. ttI 

[INe IIJ 


["AT , T l 

[INe vj 


["AT.. TTT 1 

[INe IIIJ 


rc tttI 

[o IIIJ 


["AT .. , T l 

[INe VJ 


[U IVJ 






1 n 

-lu.U-l 


1 9 81 


14 3 


A.O.OO 


1 8 71 

lO. 1 1 


94 ^1 




33 


2.8(0.2) 


1.2(0.1) 


3.9(0.2) 


2.0(0.3) 


5.3(0.2) 


2.5(0.4) 


1.6(0.2) 


8.1(0.2) 




0.022 


0.012 


0.037 


0.019 


0.051 


0.028 


0.029 


0.159 


55 


1.82(.06) 


2.2(0.3) 


< 0.5 


1.1(0.2) 


2.0(0.3) 


< 0.4 








0.062 


0.17 


< 0.01 


0.036 


0.068 


< 0.01 






172 


< 0.1 


< 0.2 


< 0.3 


< 0.6 


< 0.6 


< 0.8 


... 


... 




< 0.08 


< 0.1 


< 0.08 


< 0.2 


< 0.2 


< 0.3 






220.1 


< 0.2 


0.5(0.2) 


< 0.2 


< 0.5 


0.35(.09) 


< 0.2 








< 0.06 


0.21 


< 0.05 


< 0.1 


0.11 


< 0.05 






234 


1.7(0.1) 


3.2(0.3) 


0.8(0.2) 


3.3(0.7) 


8.2(0.7) 


< 3. 


3.9(0.7) 


7.5(1.0) 




0.0034 


0.0077 


0.0022 


0.010 


0.027 


< 0.01 


0.026 


0.056 


244.1 


0.7(0.2) 


0.9(0.2) 


1.4(0.2) 


0.6(0.4) 


0.3(0.2) 


0.68(.06) 








0.033 


0.041 


0.067 


0.033 


0.02 


0.043 






263.1 


< 0.6 


< 0.8 


0.6(0.2) 


< 0.6 


0.6(0.2) 










< 0.1 


< 0.3 


0.19 


< 0.2 


0.23 








265 


0.6(0.2) 


1.6(0.4) 


< 0.8 


< 0.6 


< 0.5 










0.012 


0.040 


< 0.02 


< 0.02 


< 0.02 








268.1 


0.49(.08) 


< 0.4 


< 0.4 


< 0.8 


< 0.2 










0.15 


< 0.1 


< 0.1 


< 0.4 


< 0.04 








280 


< 0.3 


< 0.5 


0.3(0.1) 


< 0.4 


0.9(0.2) 










< 0.01 


< 0.02 


0.014 


< 0.02 


0.050 








330 


0.4(0.2) 


0.64(.08) 


< 0.5 


0.4(0.1) 


0.7(0.2) 


0.19(.08) 








0.03 


0.076 


< 0.05 


0.04 


0.09 


0.02 






381 


1.3(0.1) 


0.62(.07) 


0.6(0.2) 


0.7(0.3) 


1.2(.3) 


1.49(.08) 


0.56(.08) 


3.9(0.2) 




0.019 


0.012 


0.011 


0.014 


0.025 


0.040 


0.026 


0.196 


433 


3.1(0.3) 


2.2(0.2) 


1.9(0.3) 


2.7(0.3) 


5.2(0.3) 


1.2(0.8) 


2.3(0.5) 


7.9(0.4) 




0.014 


0.023 


0.012 


0.020 


0.042 


0.013 


0.028 


0.105 


452 


< 0.5 


0.7(0.1) 


2.11(.07) 


< 0.3 


1.9(0.2) 


1.4(0.4) 


< 0.3 


1.3(0.2) 




< .008 


0.012 


0.034 


< .005 


0.032 


0.026 


< 0.01 


0.057 
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a Notes. For each source, emission line fluxes (10~ 14 erg s _1 cm -2 ) or 2a upper limits are 
on the first line and rest equivalent widths (/im) are on the second line. Emission line rest 
wavelengths (/im) are at the top of each column. 
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Fig. 1. — Comparison of 15 /im (rest) mid-IR power to 178 MHz (rest) radio power. The 
dotted lines indicate mid-IR to radio ratios of 1, 10, 100, and 1000. Quasars and BLRGs 
(asterisks) and mid-IR luminous NLRGs (diamonds) have mid-IR luminosities i/L„ (15/im) > 
8 x 10 43 erg s~ 4 (dot-dash line). Mid-IR weak NLRGs (open triangles) have uL u (15/j,m) < 
8 x 10 43 erg s -1 . Filled triangles denote mid-IR weak NLRGs with 3a upper detection limits 
at 15 /im. 
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Fig. 2. — Spitzer spectra of mid-IR luminous narrow-line radio galaxies (z = — 0.2), 
ordered by 9.7 fim silicate trough optical depth. The continuum fits used to measure the 
silicate troughs are shown as dashed lines. The peak at f 2-20 //m is characteristic of thermal 
dust emission. Several high-ionization forbidden emission lines, H 2 pure rotational lines, and 
the 11.3 /iin PAH feature are labeled. 
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Fig. 3. — Spitzer spectra of mid-IR luminous NLRGs (z = 0.2 — 0.7). Some data are rebinned 
(diamonds) to improve S/N. Silicate absorption troughs (below dashed lines) and forbidden 
emission lines are present in some sources. 
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Fig. 4. — Spitzer spectra of mid-IR luminous NLRGs (z = 0.7—1.0). Some data are rebinned 
(diamonds) to improve S/N. Some of the large downward spikes and apparently missing data 
for 3C 263.1 and 268.1 are artifacts of plotting noisy data on a logarithmic scale. Silicate 
absorption troughs (below dashed lines) and forbidden emission lines are visible in some 
sources. 
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Fig. 5. — Spectral energy distributions of mid-IR-luminous NLRGs with near-IR photometry 
(triangles), ordered by 15 fim luminosity. Spitzer data (solid curves) reveal excess mid-IR 
emission relative to the near-IR stellar component. Blackbody fits (dashed curves) to the 
low and high-frequency ends of the Spitzer spectra indicate dust temperatures of 210-660 K. 
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Fig. 6. — (a) Mid-IR luminosity vs. 7-15 //m spectral index for mid-IR luminous NLRGs 
(diamonds), mid-IR weak NLRGs (triangles), and quasars or BLRGs (asterisks). Mid-IR 
weak NLRGs with a < 1.0 (e.g., 3C 310) have significant starlight from the host galaxy at 
7 /im. The reddening vector corresponds to 100 magnitudes of optical (V band) extinction, 
(b) Silicate 9.7 /im apparent optical depth vs. 7-15 /im spectral index for mid-IR luminous 
NLRGs with well-measured troughs. The two reddening vectors both correspond to 10 
magnitudes of optical (V band) extinction, for two different Galactic-type dust extinction 
curves (see text). 
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Fig. 7. — High-ionization emission line/ [Ne il] A12.81 /im ratios vs. mid-IR luminosity 
vL v {lb /im, rest) in units of 10 43 erg s _1 , for mid-IR luminous NLRGs. (a) [O iv] A25.89 
/im, (b) [Ne Vl] A7.65 //m, (c) [Ne v] A14.3 /im (open diamonds), and [Ne v] A24.31 /im 
(filled diamonds), (d) [Ne ill] A15.55 /im, (e) [S iv] A10.51 /im, and (f) [S ill] A18.71 /im. 
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Fig. 8. — Mid-IR (15/xm, rest) luminosity vs. 5 GHz (observed) radio core luminosity, both 
normalized by 178 MHz (rest) radio lobe luminosity. The various symbols represent mid- 
IR weak NLRGs (triangles, filled = upper limit), mid-IR luminous NLRGs (diamonds), and 
quasars (asterisks). Upper limits at 5 GHz are indicated by arrows. The most core-dominant 
sources, with vL v {h GHz)/z/L^(178 MHz)> 10 2 are all quasars or BLRGs. 



